Composition and concentration anomalies for structure and dynamics of 

Gaussian-core mixtures 
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We report molecular dynamics simulation results for two-component fluid mixtures of Gaussian- 
core particles, focusing on how tracer diffusivities and static pair correlations depend on temperature, 
particle concentration, and composition. At low particle concentrations, these systems behave like 
simple atomic mixtures. However, for intermediate concentrations, the single-particle dynamics of 
the two species largely decouple, giving rise to the following anomalous trends. Increasing either the 
concentration of the fluid (at fixed composition) or the mole fraction of the larger particles (at fixed 
particle concentration) enhances the tracer diffusivity of the larger particles, but decreases that of 
the smaller particles. In fact, at sufficiently high particle concentrations, the larger particles exhibit 
higher mobility than the smaller particles. Each of these dynamic behaviors is accompanied by a 
corresponding structural trend that characterizes how either concentration or composition affects the 
strength of the static pair correlations. Specifically, the dynamic trends observed here are consistent 
with a single empirical scaling law that relates an appropriately normalized tracer diffusivity to its 
pair-correlation contribution to the excess entropy. 



Fluids of identical particles interacting via the 
Gaussian-core (GC) pair potential have been the sub- 
ject of many recent investigationsi^i^i^i^i^iiiSii^ Contin- 
ued interest in this model system, introduced by Still- 
inger in 1976,^" can be attributed in part to the fact 
that the GC potential is a simple and computationally 
tractable idealization of the soft, effective interparticle 
repulsions that can exist between large molecular species 
(e.g., star polymers) or self-assembled structures (e.g., 
micelles) in solutionj^i The GC fluid is also a compelling 
model to study because it exhibits several unusual phys- 
ical properties that are typically associated with molec- 
ular or complex-fluid systems with more complicated in- 
teractions. For example, at low temperature, the GC 
fluid displays a re-entrant freezing transition^i^iiii^ neg- 
ative thermal expansivity, ^'^^ and its isothermal com- 
pressibility increases upon isobaric coolingj^ Although 
the structural and dynamic properties of the GC fluid 
are qualitatively similar to those of simpler fluids at low 
particle concentrations, they become anomalous at suf- 
flciently high particle concentrations. For example, the 
single-particle dynamics, quantifled by, e.g., the self diffu- 
sivity, become progressively faster upon increasing parti- 
cle concentration (diffusivity anomaly) The static 
pair correlationSji^'^'^'^'^'i^ quantifled by, e.g., the two- 
body excess entropy s*-^-*, also weaken upon increasing 
particle concentration (structural anomaly) 

The differences between the structural behavior of the 
GC fluid at low versus high particle concentration can 
be qualitatively understood by considering the Gaussian 
form of the repulsion. At low concentration and low tem- 
perature, the average interparticle separation is larger 
than the range of the interaction. Thus, the part of the 
GC potential that the particles typically sample when 



they "collide" is steeply repulsive. Under these condi- 
tions, small increases in concentration lead to the build 
up of short-range static correlations (i.e., coordination 
shell structure) , similar to what occurs in the hard-sphere 
fluid i^i^i'^i^'^i'SiiSi'iS However, at sufficiently high concentra- 
tion, the bounded form of the GC potential allows the 
average interparticle separation to become much smaller 
than the range of the interaction. As a result, parti- 
cles are effectively penetrable and constantly experience 
soft repulsive forces from many neighbors. These forces 
largely cancel one another, creating a "mean field" . Fur- 
ther increasing the concentration only makes this effect 
more pronounced, paradoxically driving the high-density 
system toward an ideal-gas-like structurci^i^i^i^iii^iSiii^ 

Less is understood about the microscopic origins of 
the anomalous relationship between diffusivity and par- 
ticle concentration, although the results of recent in- 
vestigations indicate that the unusual dynamical trends 
are closely linked to the aforementioned structural 
anomaliesi^ In particular, the equilibrium GC fluid ex- 
hibits a semi-quantitative scaling relation^ between self 
diffusivity D and the two-body excess entropy s*-^) . In- 
terestingly, this relationship is "normal" in the sense that 
it is similar to that observed for a wide variety of simpler 
fluids that do not exhibit either structural or dynamic 
anomalies. Stated differently, the diffusivity 
anomaly of the equilibrium GC fluid disappears when 
one plots D versus s*^^^ instead of particle concentration)^ 
Similar trends have also recently been reported for other 
equilibrium fluids with dynamic and structural anoma- 
lies, e.g., models with water- like interaction a^^i^^'^^ or 
colloid-like, short-range attractionsi^i^ 

In this paper, we further explore the relationship be- 
tween structure and dynamics in simple models for com- 
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plex fluids by studying, via molecular simulation, binary 
mixtures of GC particles. The fluid phase behavior of 
these systems has already been studied extensivelyi^ 
However, here we present, to our knowledge, the first 
investigation of the relationships between the static pair 
correlations of the fluid and the tracer diffusivities of the 
two components. 

Specifically, we study the following questions about 
how these quantities depend on particle concentration 
and mixture composition. Are the trends in the tracer 
diffusivities of the two components of the GC mixture 
closely coupled? Do they scale in a simple way with a sin- 
gle measure of the overall strength of the pair correlations 
(e.g., s^^^^)? Or, alternatively, is there a significant decou- 
pling of the single-particle dynamics of the two species? If 
this latter scenario holds, do the resulting trends in tracer 
diffusivities track decoupled, species-dependent measures 
of static structure? Finally, what are the implications of 
the answers to the above for the compositional depen- 
dencies of structural order and tracer diffusivity at low 
versus high particle concentration? 

To address these questions, we use molecular dynam- 
ics simulations to investigate equilibrium two-component 
fluid mixtures of particles that interact via pair potentials 



of the GC form, 4>ij{r) = exp[— (r/tTi 



Here, r is 



the interparticle separation, and the parameters eij and 
Uij characterize the energy scale and range of the interac- 
tions, respectively, between particles of type i and j with 
J, j G {A, B}. Since we want to understand the behavior 
of uniform binary fluids, we assign numerical values to 
the parameters that favor mixing. Speciflcally, we adopt 
a set^ introduced earlier {ubb — 0-665<7aa', o^ab = 
(0.5[cr^^ + cr|B])°-^ eAA = esB] ^AB = 0.9446^^), in 
which the Uij were chosen to mimic binary mixtures of 
self-avoiding polymers in solution)^ We truncate all pair 
potentials at an interparticle separation of •i.2aAA, and 
treat the particles of the two species as having equal 
masses {mA — tub)- 

We carry out the simulations in the microcanonical 
ensemble, numerically integrating Newton's equations of 
motion with the velocity- Verlet scheme^ using a time 
step of Qm^mA(T\JeAA- We use N = 3000 GC par- 
ticles and a periodically replicated simulation cell, the 
volume V of which is chosen to realize speciflc values of 
reduced total concentrations (i.e., particle densities) in 
the range 0.05 < pct\a ^ 1; where p = N/V. We inves- 
tigate mixtures over a wide range of composition (0.1 < 
XA < 0.9, where xa is the mole fraction of species A) 
and reduced temperature (0.05 < kBT/tAA < 0.40). To 
characterize the single- particle dynamics of species i, we 
compute its tracer diffusivity Di by fitting the long-time 
(t oo) behavior of its average mean-squared displace- 
ment (Ar^^) to the Einstein formula, = (Ari^)/6ii^ 

We compute the two-body excess entropy s^^^ directly 
from the partial radial distribution functions gij (r) of the 
fluid using the expressiou j^^i^^ 




FIG. 1: (a)Tracer diffusivity Di and (b) structural order met- 
ric -s^f \ with i G {A,B}, versus concentration pctaa for the 
binary Gaussian-core fluid mixture discussed in the text. The 
collective structural order metric -s^^^ is also included in (b). 



The temperature is ksT/eAA 
XA = 0.5. 



where s|^^ is given by 



0.1 and the mole fraction is 



(2) 

fee ^ 



[dz] (r) In g^J (r) - (r) + l]dr (2) 



(1) 



Both — s^^^ and — s.p^ are non- negative and can be viewed 
as translational structural order metricsi^ The former 
characterizes the overall strength of the pair correlations 
in the fluid, while the latter quantifies the amount of 
pair structuring surrounding particles of type i. 

The first issue that we investigate using our simula- 
tion data is how closely the single-particle dynamics of 
the two species are coupled. Figure [TJa) shows how the 
computed tracer diffusivities. Da and Db, depend on 
density pcr\j^ for an equimolar {xa = 0.5) mixture at 
a temperature of k^T/eAA = 0.1. As can be seen. Da 
follows the same type of non-monotonic trend observed 
for the self diffusivity of the single-component GC fluid, ^ 
displaying an anomalous dependency on particle concen- 
tration [{dDA/dp)T,XA > 0] for densities greater than 
P(7aa ~ 0-4- On the other hand, Db shows behavior 
consistent with that of simple fluids, monotonically de- 
creasing with per A A over the density range examined here. 
The fact that Da and Db decouple in this way gives rise 
to a dynamic crossover density, above which the larger 
A particles exhibit higher mobility than the smaller B 
particles. It also suggests that one cannot trivially corre- 
late both the Da and Db trends with a single, collective 
measure of structural order for the fluid, such as s^^^ [see, 
e.g., Fig.[T](b)], a point we examine further below. 

Does increasing particle concentration result in a corre- 
sponding decoupling of species-speciflc structural metrics 
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that, in turn, correlate with the dynamical trends of A 
and B particles? To examine this possibility, we first plot 
in Fig.[TJb) the density dependencies of —s\ and — 
As can be seen, there is indeed a structural decoupling. 
While both 



-s^'' and 



increase with density at low 
values of /ocr^^ , the increase in ordering is initially more 
pronounced for the structure surrounding the A particles. 
This is to be expected because A particles have larger ef- 
fective contact diameters with their neighbors, and thus 
respond by building up stronger static pair correlations 
at low density. However, the large size of the A particles, 
coupled with the bounded nature of the GC interaction, 
means that A particles are also forced into more inter- 
particle overlaps than the B particles as the density is 
increased. This ultimately leads to a weakening of the 
static correlations (i.e., coordination shell structure) of 

(2) 

the A particles, and hence a maximum in —s\ at an in- 
termediate value of P(t\a- As should be expected based 
on the smaller size of the B particles, a maximum in 

— (and a corresponding minimum in Db) can also 
occur at significantly higher densities, if phase separa- 
tion of the mixture does not occur first. ^^Interestingly, 
similar to what is observed for the tracer diffusivities, 

(2) 

one of the consequences of the s,- decoupling described 
above is the presence of a structural crossover density 
above which the smaller B particles exhibit stronger pair 
correlations than the larger A particles. 

The data in Fig. [T](a) and (b) suggests a negative corre- 



,(2) 



(and also between Db and 



lation between Da and 

~s^^), in which the structural crossover (s^'' = s'^b) 
curs at approximately the same density as the dynamic 
crossover {Da = Db)- In fact, although we focus on 
one particular pairing of composition and temperature 
in Fig. [H these trends are exhibited by this system for 
a wide range of compositions and temperatures (see the 
supplementary Fig. 2]) 

To quantitatively examine the correlation between 
single-particle dynamics and structure, we studied possi- 
ble generalizations of a scaling relationship that describes 
the behavior of the single-component GC fluid. In the 
case of the single-component fluid, the so-called Rosen- 
feld scaledi^ self diffusivity D^ = Dp^/^{kBT /rrii)-^/^ 
is approximately a single- valued function of — s'^' across 
a wide range of temperature and densityi^ Figure [H[a), 
however, shows that a nai've extension of this result 
for mixtures, Df = Dip^/^{k^T /mi)~^/'^ versus —s^'^\ 
does not adequately collapse the data for either of the 
two species. This should not be particularly surprising, 
given the dynamical and structural decouplings shown in 
Fig. [Ha) and (b). 

On the other hand. Figure HUb) examines a species- 
specific extension of the single-component scaling law, 
Df versus ~s'f'\ Interestingly, not only does this gener- 
alization collapse the temperature, density, and compo- 
sitional dependencies of tracer diffusivity for each indi- 
vidual particle type, but the behaviors of the two species 
are, to a good approximation, accounted for by the math- 
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versus (a) two body excess entropy 



s'^-* and (b) its contribution from structuring around type i 
particles -s^^ with i £ {A,B}, for the binary Gaussian-core 
fluid mixture discussed in the text. Shown are mole frac- 
tions xa ~ 0.1 (black), 0.3 (red), 0.5 (orange), 0.7 (blue), 
and 0.9 (green) and temperatures kBT/eAA = 0.05 (o), 0.1 
(□), 0.2 (o) and 0.4 (v). Filled and open shapes represent 
A and B particles, respectively. The dashed line indicates 
a least-square fit to a power law relationship for the single- 



component GC fiuid?, D 



■ 0.208s(2)" 




FIG. 3: Tracer diffusivity Di for (a) A particles and (b) B 
particles and structural order metric -sp' for (c) A particles 
and (d) B particles for the binary Gaussian-core mixture dis- 
cussed in the text. All are plotted versus density p<j\a- Data 
shown is for temperature k^T/eAA ~ 0.2 and mole fractions 
Xa ~ 0.1, 0.3, 0.5, 0.7 and 0.9. Arrows indicate increasing xa- 
The dashed lines separate the approximate low, intermediate 
and high density ranges described in the text. 
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ematical form of the scaling law for the single-component 
GC fluid. 

We now explore what this correlation implies about 
how mixture composition affects the tracer diffusivities 
[Figure[3fa),(b)] and species-specific pair-correlation con- 
tributions to excess entropy [Figure [3l^c),(d)]. As should 
be expected, at low values of density (pcr^^ < 0.2), the 
response of the system to changes in composition is nor- 
mal, i.e., qualitatively similar to that of simple atomic or 
hard-sphere-like mixtures where interparticle "collisions" 
dominate. Under these conditions, increasing the mole 
fraction of the larger A particles effectively increases the 
"packing fraction" of the fluid, which in turn decreases 
the mobility and increases the local structural order sur- 
rounding both types of particles. 

At high values of density, on the other hand, the fact 
that the bounded GC potential allows for significant in- 
terparticle overlaps changes the physics. Since increas- 
ing the mole fraction of the larger particles (at con- 
stant density) here increases the number of overlaps and 
nudges the system toward the mean-field fluid, one ex- 
pects anomalous behavior in dynamics and structure, i.e. 
a corresponding increase in Di and decrease in —sf'^ for 
both species. Indeed, Figure [3] shows that these anoma- 



lous compositional trends for dynamics and structure do 
occur for pcr^^ > 0.4. 

In closing, we note a final manifestation of the de- 
coupled behavior for the two species of this GC mix- 
ture. Specifically, the densities at which the composi- 
tional trends for structure and dynamics transition from 
normal to anomalous are significantly different for the 
two species, with the larger species logically becoming 
anomalous at a lower overall density. The implication 
is that there is a fairly wide range of intermediate fluid 
densities (approximately 0.2 < p(j\j^ < 0.4) for which 
the structure and dynamics of the A particles behave 
anomalously, while those of the B particles behave nor- 
mally, with respect to changes in composition. 
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FIG. 4; (a-d)Tracer diffusivity Di and (e-h) structural order metric -s] , with i G {A,B}, versus concentration pcjaa for the 
binary Gaussian-core fluid mixture discussed in the text. The columns are systems at different temperatures and mole fractions, 
which are listed above the respective column. 



